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ABSTRACT 

This paper demonstrates analytically how errors in tracking 
station time and frequency synchronization, as well as errors in 
station location, influence the accuracy of trajectory determina- 
tion. Two systems, whose data are presently used for most of the 
more accurate orbit determination schemes, are described: a 

range and range rate system, and a radar system. All other sys- 
tems are actually combinations of these two basic systems. 

Rather simple analytical expressions are derived relating 
measuring errors with those of time and frequency synchroniza- 
tion, as well as ground tracking station location. For the range 
and range rate system, the error in range rate Sr is used as a 
yardstick, and all other quantities are derived from it. For the 
radar system (not measuring range rate directly), the total local 
position vector, as measured by such a system is used as the 
yardstick to determine the necessary time synchronization and 
station location accuracy. Numerical examples are presented 
which are intended to show what accuracies of the mentioned 
quantities are needed for a good ground tracking system to make 
data most useful in determining accurate orbits and space 
trajectories. 
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SPACE TRAJECTORIES AND ERRORS IN TIME, 
FREQUENCY, AND TRACKING STATION LOCATION* 


by 

F. O. Vonbun 

Goddard Space Flight Center 


SUMMARY 

General 

Frequency— or more definitively, phase change over a period of time— is among the few physi- 
cal quantities that can be measured with extreme accuracy. If, therefore, such measurements can 
be used to determine orbits, one would expect the determinations of such orbits to be very accurate. 
This is in fact the case, and Doppler- frequency measurements are so used on an almost daily basis, 
particularly for computing trajectories in deep space. In brief, the Doppler frequency component 
in an electromagnetic signal from a spacecraft is often used as a yardstick. 

One could use, even more appropriately, an error quantity at a particular target: for instance, 
the closest approach to the moon or a planet should be, say, 100 meters or 10,000 meters. How- 
ever, using an error quantity in this manner may necessitate system requirements that cannot be 
met because of physical limitations. These limitations include systems noise, bias errors, and 
tropospheric and ionospheric propagation that will disturb the transmitted and received electro- 
magnetic signals. 

For this is why it has been decided here to limit the range rate. In any event, the final results 
should be almost the same. Since this yardstick is an extremely accurate one (a^/V^ 10“ 10 to 10~ 12 
per day), all other quantities involved must be appropriately accurate. For example, for the Dop- 
pler measurement to be exploited to its fullest, the sender frequency must be very accurate, and 
the time synchronization between tracking stations, as well as their location on the surface of the 
earth, must be determined to a high degree of certainty. This paper outlines how much and in 
what ways these quantities and their errors are approximately related to each other. Coherent 
two- or three-way, phase- locked Doppler systems as well as radar systems are also discussed. 

Establishing reasonable design limits cannot be mathematically rigorous. Human experience 
and judgment are important factors. Nevertheless, it is felt that the numerical values arrived at 
will give the reader a good picture of the problem. 

* Published also as NASA Technical Memorandum X-55798, 1966. 
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Relative-Frequency Stability 

As can be noted later on (from Equation 9) one needs a relative frequency stability of 



2 x 10” 9 , 


where the first term refers to the measured Doppler shift during the measuring time t m , the sec- 
ond term refers to the relative stability of the local standard oscillator over the time the signal 
travels from the station to the spacecraft and back, Sr =0.01 cm/sec (the assumed yardstick), 
and r is assumed to equal 50 km/sec for an average maximum interplanetary range rate. In addi- 
tion, the value (Sc/c) must be considered as an unknown in the orbit determination process in 
order that the above values of Sr 3 etc., can be used. 

Station Time Synchronization 

Equations 12 and 13 (see further on) give St = 0.1 to 0.2 msec for the earth- moon space using 
Si- = 0.1 cm/sec and full earth acceleration p = 1000 cm/sec 2 (a pessimistic value proper for 
near the earth). 

Similar St =0.2 msec is the needed time synchronization error for a radar depending on the 
assumed angular bias error of the radar considered. Equations 16, 17, and 17a are applicable for 
this case. 

For interplanetary missions where p is small (in this case p = 0.6 cm/sec 2 is the accelera- 
tion due to the sun for 1 AU), even with Sr = 0.01 cm/sec (compared to 0.1 cm/sec for near-earth 
tion due to the sun for 1 AU), even with Sr =0.01 cm/sec (compared to 0.1 cm/sec for lunar 
space), a value of St =2 msec is adequate as given by Equation 13. Only in the special case 
of a planetary approach or a planetary fly-by, timing errors of St = 10 to 20p sec seem to be 
the orbit for time synchronization, can fulfill this requirement without additional equipment or 
methods. 


Station Location Errors 

To make full use of a Doppler error as small as 0.01 cm/sec requires a station location ac- 
curacy of: 


| SR | - 1.5 to 5.5 meters , 


depending on the use of Equation 26 and 27. This range of accuracy is something to aim at. 


2 



In the case of a radar system, a system error range of 


| 6R\ - 3. 5 to 11.0 meters 

will be necessary for future systems according to Equation 30 and its associated assumptions. 


Frequency Synchronization 

For a two-way Doppler system, similar values hold— as stated under Relative- Frequency 
Stability above. In this case one speaks not of frequency synchronization as such but, rather, of 
frequency stability during the travel time of the electromagnetic signal. (This is, of course, dif- 
ferent for the three-way Doppler system, as applied to lunar and planetary missions.) 

Equation 19, shows that to make full use of the Doppler (Sr = 0.01 cm/sec) requires a fre- 
quency synchronization between stations 1 and 2 of 


61.12 




0 



Or, in brief, field- worthy hydrogen masers must be developed and used. 

Note that the errors considered— namely, the errors in frequency synchronization, tracking 
station time synchronization, and location of the tracking stations on earth— are to be considered 
as bias errors. This is obvious, since all these quantities stay approximately constant while track- 
ing measurements are taken and the orbits are determined. 


FREQUENCY STABILITY 

One of the first questions to be answered in 
designing a coherent Doppler tracking system 
is: what frequency stability of the transmitter 

is required? 

Figure 1 shows the principle (References 1 

through 5) used to detect the Doppler frequency 

shift for a two-way Doppler system. 
u 2 

The multiplied frequency of a stable 
source is transmitted to the spacecraft, re- 
ceived there, translated to ensure no inter- 
ference, retransmitted to the same ground re- 
ceiving system together with the translation 
factor k, and finally fed into a mixer (Reference 
5). Here, the difference frequency between the 


SPACECRAFT 



Figure 1 —Two-way Doppler system. 
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signal sent to the spacecraft and that received back on the ground is extracted as the Doppler fre- 
quency v D 2 ; this is trustworthy, provided that the oscillator maintains constant frequency during 
the signal travel time, 2 r. The significance of "constant" will be shown later. 


The Doppler Shift 

The special relativistic equation for the Doppler shift reads (References 1 through 4): 


1 _ 

1 - ^-(v • r 




( 1 ) 


where v is the relative velocity between spacecraft and receiving station (v= |v|); c the speed 
(scalar) of light; v the frequency at the spacecraft; v 0 the transmitted frequency, and r° the unit 
position vector from the station to the spacecraft in the moving frame of reference. Further, 


r - ( v * r°) 


( 2 ) 


is the spacecraft range rate with respect to the tracking station. A similar equation holds for the 
return path, as shown in References 1 through 4. The ratio of the received frequency v 1 and the 
transmitted frequency is shown to equal 


The Doppler frequency for the so-called two-way mode v D is known simply as: 

"d 2 = ■ ( 4 ) 

The factor k can be assumed to be 1 since this value is compensated for on the ground, as shown 
in Reference 5. As mentioned, the frequency translation (k) was performed only to eliminate in- 
terference at the spacecraft. 

Equations 3 and 4 give for the range rate 


* c 


1 + c 


(3) 


dr 

dt z 



1 

1 + ( v d/ 2 v o) 


(5) 


Equation 5 relates Doppler frequency ^ D , transmitter frequency v 09 speed of light c, range rate r 
between spacecraft and transmitting stations in exact form as given by the special theory of rela- 
tivity. (No gravitational terms are considered since their effects are relatively small.) 
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The Relative Errors 


Equation 5 will now be used to answer the question of frequency stability using (dr/r)* ** as a 
"yardstick. M Since only small variations are to be considered and v D /v 0 «< 1, Equation 5 can be 
simplified for the following variational analysis, that is: 


Varying Equation 6 and adding the results in the Gaussian sense leads to 




(7) 


This equation shows that the relative error Sr/r in the determination of r depends on the relative 
error in the speed of light Sc/c, the relative error in the Doppler frequency §^ D /^ D , and the trans- 
mitter frequency . Decreasing Sr or b r/r means reducing all the other three quantities. 

First, it is necessary to improve the relative error (Sc/c) of the speed of light (at present 
(Sc/c) = 3 x 10” 7 (References 6 through 9), then measure the relative Doppler frequency (s ^ D /^ D ) 
to at least the same accuracy. This means that: a certain measuring time t m , and measuring 
time error St m must be secured to make a frequency measurement to a given accuracy (Reference 
10), that is: 


Sz'n 


St 

IT 

__ 


2r 

c 


St 

IT 


(8) 


assuming that the quantities contributing to noise (noise modulation index, signal to noise ratio) 
are small enough compared with (St m/ /t m ), which is true for the systems and methods considered 
here. 

For example, a carrier frequency = 2000 Me, a measuring time t m = 1 sec, and an error 
in the measuring time St m = 10" 6 sec give a value of 

Sl/ D 2 2 

— - “3 * 10 4 ' 10” 6 * - ■y 10“ 10 * (for r = 10 km/s) . 

In addition, of course, the transmitter standard oscillator must be "constant” during the travel 
time 2r of the electromagnetic wave from the transmitter to the spacecraft. A shift during this 


*During the course of this paper normalized errors such as S^c, S^r. . . . etc., are often used for simplicity of the mathematical 
expressions. 

**Sign of Doppler frequency is neglected from now on since a "negative/* frequency has no real meaning. 
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TIME FROM EPOCH (hours) TIME OF RCA 
Figure 2— Position errors of a lunar transfer trajectory. 



TIME FROM EPOCH (hours) f 


TIME 
OF RCA 


Figure 3 — Velocity error of a lunar transfer trajectory. 


time cannot be detected by the mixing process 
shown in Figure 1 and described in Reference 5. 
In brief, one needs: 



Assume, as an example, that all three quantities 
equal 10~ 6 (a conservative value) than Sr = V3.10“ 6 r. 
For flights between the earth and the moon, r ar 
= 3 km/sec = 3.10 s cm/sec and thus Sr = 1.0 
cm/sec which is at present adequate for pur- 
poses of trajectory determination for the Apollo 
project (References 7, 8, and 11 through 16). 

(See also Figures 2 and 3.) Figure 2 depicts the 
instantaneous position error (propagated to the 
moon) of a lunar transfer using a bias range- 
rate error of 1 cm/sec plus a 1 cm/sec noise 
component superimposed. Station errors as 
given in References 7 and 8 are also included to 
show realistic values as used for these Apollo 
Navigation studies. Figure 3 shows the same 
for the spacecraft velocity. 

For interplanetary flights one must decrease 
the errors in range rate, in order to get a more 
reasonable trajectory. Figures 4 and 5 show the 
position and velocity errors of a Jupiter probe 
as an example when tracked by a two-way Dop- 
pler system using three ground stations approx- 
imately 120 degrees separated in longitude. The 
errors of twenty days of tracking are shown; 
these errors are then propagated to 260, 458, 
and 500 days respectively (planets intercept) 
along the trajectory. This is in order to be able 


to find if it is reasonable to make a midcourse maneuver after, say, 20 days of tracking based 


upon the range- rate accuracy requirements discussed here. 


Since at present both quantities— (8 v d /i/ d ) and (Sv 0 /v 0 )— can be determined far better than 
(Sc/c) one must improve on the latter to improve (Sr/r) as shown by Equations 7 and 7a. This 
should possibly be done by an orbit- independent method in a laboratory, since c and thus Sc/c can 
be calculated (considered as an additional unknown in the orbit computation) together with the orbit 
parameters station locations, etc. The problem involved here is that all our constants are based 
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upon the velocity of light within 3 x 10" 7 as 
quoted before (Reference 17). The velocity of 
light is a primary physical constant and oc- 
cupies a key position in astronomy (References 
17 and 18). The velocity of light is involved via 
the ’Tight time” (time the light needs to trans- 
verse 1 AU, the heliocentric distance of the 
center of mass of the earth + moon system) with 
the solar parallax, the mean earth radius, the 
eccentricity of the earth orbit, the Gaussian 
gravitational constant, etc. This means that all 
these quantities must be changed accordingly to 
adjust our physical constants in the universe in 
terms of distance, distance variation, time, etc. 
(Reference 17). 
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For interplanetary flights (References 19 
through 21), values of Si =0.02 cm/sec have 
been obtained over 1-minute of measuring time. 
Thus, using 0.01 cm/sec as a standard, one 
would need for r = 5. 10 6 cm/sec the following 
accuracies (extreme values): 




Figure 4— Position errors of a Jupiter trajectory. 
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Figure 5— Velocity errors of a Jupiter trajectory. 


On the assumption of a good signal- to- noise 
ratio of 10 db or more, and Doppler measuring 
times on the order of minutes (Reference 10), 
atmospheric disturbances seem to be the only 
limiting factors (References 19 through 22). As 
shown in Reference 22, this limit is of the order 
of 0.01 cm/sec rms when elevation angles € = 10 
degrees are used and daily tropospheric cor- 
rection terms are applied in the order of SN = ±5 where N = (n - 1) * 10 " 6 , and n is the index of re- 
fraction. The only unknown quantity is (Sc/c) to the accuracy stated, thus a better determination 
of Sc/c = 3. 10“ 7 as shown in Reference 6 is needed. Laboratory measurements of the velocity of 
light will probably not reach an accuracy of lp 10 10 or better in the near future. One way to cir- 
cumvent the errors of the known value of c is to introduce this quantity as an "unknown" into orbit 
determination (Jet Propulsion Laboratory). 
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STATION TIME SYNCHRONIZATION 


One of the next logical questions to answer is: To what accuracy should the time be synchron- 
ized between tracking stations used for orbit determination? 

A relationship between the time error St and the tracking quantities to be measured will be 
derived in the following. A CW range- rate system and a radar system are considered since their 
time synchronization requirements differ according to their different "yardsticks." (A radar does 
not "measure" r as such, for instance.) 

Synchronization Between Station Time and Orbit Time 

For orbit determination, the local or station measurements r, r , etc. (or any other measure- 
ment) are made at a certain time t. These values are then transmitted from the station to the 
computing center and there used for orbit determination. Thus, the station clock must be syn- 
chronized to the computing center "clock" to within a certain limit, say St. What this value should 
be in order to make full use of our standard Sr is shown in the following. 

Considering only small deviations within a CW range- rate system, one can write to a first- 
order accuracy for the variation r the following: 


Sr - r St , 


( 10 ) 


since 


.. _ S . 

r = St r ’ 



t - TIME OF MEASUREMENT 

%- PROPER POSITION VECTOR FOR SAY f = 0 AND NO 
SYNCHRONIZATION ERROR 


Figure 6— Radar measurement of slant range vector. 


where V is the time derivative of the range rate 
r and 8t is the time error to be determined here. 

The spacecraft inertial position vector (Fig- 
ure 6) p is 

p - R + r (11) 

and acceleration are: 

p - R + r (12) 

Borrowing terms from vector analysis, that is, 

R = (wx R) and further R - wx(SxR), one obtains 

~p - px(wxR)] + 7 (12a) 
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Using Equations 12a and 10 one obtains approximately: 


st = 



( 13 ) 


Here one has to remember two points, namely: the range rate error Sr was previously assumed to 
be 0.01 cm/sec, on the basis of deep- space mission results, and the range rate is not the only 
measurable quantity used for orbit determination for near-earth satellites except for the TRANET 
spacecraft. From experimental data such as Goddard spacecraft as well as lunar orbiter evalua- 
tion, a value of Sr = 1 cm/sec emerges. Improvements are under way, so that an assumption of 
dr = 0.1 cm/sec for earth and lunar space satellites seems reasonable for the near future. Using 
this value, namely dr = 0.1 cm/sec and p = 10 3 cm/sec 2 as the near-earth and lunar- space yard- 
stick, one obtains, using Equation 13, 


St - - 10“ 4 sec - 100 psec . (13a) 

The value I R | = 3 cm/sec 2 has been neglected, as compared with p = 1 g = 1000 cm/sec 2 . 

For interplanetary travel, where the acceleration terms are small (in the order of a few cm/ 
sec 2 ), timing errors can be appropriately large as dictated by Equation 13. 

The value given in Equation 13a should then suffice for most space tracking problems. 

The only time where more stringent requirements are to be considered is during a planetary 
approach. For such cases, the value given in Equation 13a may have to be eventually reduced to 
say 10 to 20 psec (p = 1000 cm/sec 2 , Sr m 0.01 cm/sec). Here, on the other hand, single tracking 
stations are mostly involved for approximately 8 hours a day, if necessary. The station time 
synchronization method suggested by JPL (Reference 27) can be applied, reducing the timing er- 
rors to the above-mentioned values without elaborate additional station equipment and/or other 
methods (e.g., VLF). 

A radar system for near-earth orbit calculations needs more than r information. Angular data 
(azimuth, elevation, or equivalent) and range data (radar) are used to determine orbits for manned 
flight programs (References 7, 8, 11, and 12). For these programs, r and Si will not serve as 
yardsticks; and yet time errors do enter the orbit analysis. (Their effect here is not severe, as 
will be shown.) 

Suppose (to assume the worst) that a near-earth spacecraft moves at 10 km/sec (during a lunar 
transfer orbit) and that the radar system would "see" a velocity of almost the full 10 km/sec. What 
would the time accuracy have to be to make the position measurement (r, a vector quantity) accur- 
ate to a certain predescribed number of meters (for example, 5 to 10)? 
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In this case, one can write (see Figure 6) 


r = r 0 + Vt , (14) 

where r is the measured slant range vector r at t = t (range r, azimuth a, elevation e) r 0 is the 
initial position vector at t = 0, and v is the velocity of the spacecraft relative to the radar sys- 
tem. If the stations are synchronized completely, the value r = r at t =0. This of course, can 
never be achieved. Therefore, varying Equation 14 yields 

Sr = Sr Q + Svt + vSt . (15) 

Considering all these errors one obtains, in the Gaussian sense, rearranging the terms and in- 
troducing r, a, and e: 

v 2 St 2 = Sr 2 + r 2 (Se 2 + Sa 2 ) + Sr Q 2 + r Q 2 (Se Q 2 + Sa Q 2 ) + Sv 2 t 2 


or 


St - ~ y^Sr 2 + 2r 2 (Se 2 + Sa 2 ) + Sv 2 t 2 . (16) 

Since Sr 2 and Sv 2 t 2 are always small compared to the middle term, (t probably less than a milli- 
second), one can write 


St = 


(17) 


where S^ B is the bias error in elevation and azimuth assumed to be equal for reasons of simplicity. 
Equation 17 thus gives an estimate of the time synchronization needed for a good radar system. For 
example, assume a near-earth orbit with v = 8000 m/sec, a radar slant range r = 200 km, and a 
elevation and azimuth bias 8js* b = 0.02 mrad; then St = 0.5 msec. 

For near-earth orbit determination, the angular measurements of a radar system play a major 
role, especially when only short tracking times on the order of 60 to 100 sec are available. In this 
case, one obtains 


st 






h_ 

v ’ 


(17 a) 


where Se B is the bias error in the random elevation angle, v is the orbital speed of the spacecraft, 
h is its height, and € max is the maximum angular acceleration. Equation 17a represents a worst-case 
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overhead pass. Using again as an example v =8000 m/sec, h = 200 km, Se B = 0.02 mrad, gives 
St = 0.25 msec. Thus even for sophisticated, and properly calibrated radars, a true synchroni- 
zation of 0.2 msec is adequate, since range, azimuth, and elevation measurements play a role in 
the 'needed time synchronization, as Equation 17 indicates. The angular errors push the needed 
time synchronization to a larger value. 

In Reference 26 it is shown how an error in station time synchronization influences the errors 
in the final Apollo orbits obtained. In general, under present assumptions, Apollo systems (Refer- 
ences 7 and 8), timing errors of approximately 10 msec do not appreciably influence the earth- 
parking and lunar transfer orbits. This is different for lunar orbits, where 10 msec do make the 
error ten times as great as with 1 msec. See Figure 7 (taken from Reference 26). 

A simple and elegant method for time synchronization was suggested in a JPL report (Reference 
27). The range is measured from two stations using a CW ranging system. From one of the ranges 
the "time" of one of the stations can be determined with respect to the other. The error in time 
synchronization bt - (i/ C ) s pos, where 8 pos is the position error of the spacecraft and c is the speed 
of light. Even a large position error, say 3 km (see Figure 5), results in St = (1/3) x 10" 8 x 3.10 3 
= 10" 5 = 10 Msec, a rather small time- synchronization error. 
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Figure 7— RMS velocity errors during lunar orbital phase (measurement biases updated). 
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FREQUENCY SYNCHRONIZATION 


As outlined under Frequency Stability, the frequency of a transmitting system must be constant 
only to a specified value during the travel time of the signal. 

This is not true, of course, when the trans- 
mitted signal is received at a "secondary" sta- 
tion as well as at the transmitting site. This 
technique, the three-way Doppler (a "pseudo Dop- 
pler") is applied for orbit determination of the 
lunar orbiter and will be used for the lunar- 
orbit determination during the Apollo mission 
(References 7, 8, and 11. See also Figure 8.) 

Using exactly the same approach as out- 
lined under Frequency Stability, one can write 
an equation similar to Equation 9, but 









i 





*02 



STATION 2 T 

1^02 




a < 


i/",k 


hv T 


hv 


12 
0 1 


8 v 


12 
0 2 


S i* 


= 2 x 10' 


(18) 


Figure 8— Three-way Doppler system. 


where the pseudo Doppler frequency (References 
7 and 8) at the second station is v D3 (three-way 
Doppler) and 8v 12 m (^ 0 i _v o 2 ) is the difference 
frequency (error in frequency synchronization) 
between stations 1 and 2 (see Figure 7) and r ps 
Equation 9 must be modified, since Sv 12 (un- 


is the pseudo range rate between the two stations, 
known) will be included in the three-way Doppler as can be seen by inspecting Figure 7. From 
Equation 18 one obtains for 8v 12 a value similar to Sv D3 . Obviously when ^ D2 is measured, the 

D2 , namely S^ D3 . Equation 18 can be written as 


error 8v l2 is contained in the error of 


Sf 




(19) 


where r ps represent the three-way (pseudo) Doppler. Example: 


then 


is an extremely small value. 


8 r 

ps 

= 0. 

1 cm/sec 

v o = 

2 x 

10 9 cps . 

5v 12 

2 

■ 10' 11 

u o 

” 3 
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This number is obtainable with hydrogen masers (References 23 through 25). Thus a three- 
way Doppler system (References 7 and 8) seems to be one of the first practical systems requiring 
a frequency standard of hydrogen- maser quality. A frequency difference hv 12 between these two 
standards still results in a bias error in r ps , which is more damaging for the orbit analysis (less 
accurate) than a random error. A bias error is the same, no matter how many measurements are 
made; a random error decreases with N" 1/2 , the number of measurements (References 13 through 
15). 

STATION LOCATION ERRORS 

Possible location errors, compatible with a range rate system and the "yardstick Si-,” as well 
as with a radar system, will now be discussed. These errors are considered a limit for these 
systems, based on their accuracy assumptions. 


For a CW Range-Rate System 

The aim here is to relate a station deviation SR (the variation of the station position vector 5 ) 
to the range rate deviation Sr. 


Range rate f, a scalar quantity, is simple; 
the projection of the spacecraft velocity v minus 
the station velocity R onto the unit local satel- 
lite position vector r ° (see Figure 9). In vector 
notation, this means 

r = (r° • (v-r)) = T (? • (v -5)) . (20) 

Varying Equation 20 yields the deviation Sr as a 
function of the station position error I S§| . 

Before this is done, Equation 20 will be re- 
grouped for ease of calculations and the value 

IT 3 = S = (Sx8) (21) 

borrowed from basic vector analysis intro- 
duced. The quantity 3 is the rotation vector of 
the earth (|3| = 7.3 10~ 5 sec" 1 ). 


T=V 



T - SLANT RANGE VECTOR 
7 - SATELLITE POSITION VECTOR 

Figure 9— Tracking station-satellite geometry. 
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Equation 20 can now be written as 


(23) 


[(r * v) " (r • (SxR))] 


and varying yields after some manipulation 

rSr - ^>R * (v - (oo x R)) -rr° +(rxw)| (24) 

with the assumption of Sv = 0 and = 0. In order to "see” the influence of the station location 
error SR onto the range rate f, one has to assume a perfect orbit, that is Sv = 0. In addition, the 
variation Sw of the earth rotation is certainly zero during the short time considered. 

Regrouping Equation 24 yields 


or 



Sr - “ uSR cos x rl + aSR • sin<5£cos ^ 2 


Since the angles 4 \ 2 are not known, their statistical average value 



1/2 


cos 2 x dx 


1 

Y2 


(25) 

(26) 


is used. For the velocity |u| T |v-r°f! the value (v 2 -r 2 ) is used as an approximation. Thus, 
Equation 26 finally reads 


Sr 



u 2 sin 2 


(27) 


where <t> is the angle between the satellite position vector r and the vector of the earth T s rotational 
axis a-. 

Equation 27 thus shows the relationship as mentioned at the beginning of this section. This 
equation gives as aimed at a rather simple relationship in order to obtain a ’’feel" for the situa- 
tion and to estimate the error of r that is sf. 

Example: (1) Deep space mission (Reference 19) 

dr ~ 0.01 cm/s 


r — very large or 
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SR 


fi • Sr_ _ 1 . 4x 10' 2 

sin 4>’ <» 0. 63x 7.3x 10' 5 

= 300 cm = 3 meters . 


(2) Near earth mission (Reference 8) 


Sr = 3 cm/s = 0.03 m/s 

v - 8 km/s 

r = 6 km/s 

r - 500 km 

then o 2 sin 2 <£ can be neglected in Equation 27 so that 


S r 


SR & ~ r 2 

72 


or 


also 


r 



- 500 

1.4 3 10 2 ^=r - 4 meters . 

^64 - 36 

Please note that considerable different range rate errors (Sf = 0.01 cm/s, 3 cm/s) were used 
based upon past and present experience for these particular missions. 

Nevertheless, it shows that station location is important if one intends to make ’'full use” of 
the range rate information. 

For a Radar System 

Suppose that without measuring r as such (another criterion) the position error must be in- 
troduced in a radar system. In general, if the station position error SR is smaller than the radar 
error Sf , no improvements can be obtained any longer by the radar system. That is 


SR = Sf 


(28) 
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or 


SR = /sr 2 + r 2 (S« 2 + Sa 2 ) 


(29) 


Using the same reason (neglecting Sr 2 ) as for Equation 17, one obtains for the station error 
(square root of the number of the components) 


SR 


ft 





(30) 


where 8*. = Sa = Se, N the number of measurements used to Construct" one radar point (a vector, 
?(r, a, e)) and S\ B is the angular bias error. 

Example: Let N = 10, r = 300 km, = 0.02 mrad, = 0, and 0.02 mrad. Then we obtain: 

SR = 11 meters 


and 


SR - 3.5 meters , 

depending on the bias error S^ B . These figures are believed to give a reasonable range of station 
accuracies. 
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